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Abstract. Accurate numerical solution of the five-body Schrodinger equation is effected via vari¬ 
ational Monte Carlo in a correlated quark model. The spectrum is assumed to exhibit a narrow 
resonance with strangeness S = + 1. A fully antisymmetrized and pair-correlated five-quark wave 
function is obtained for the assumed non-relativistic Hamiltonian which has spin, isospin, and color 
dependent pair interactions and many-body confining terms which are fixed by the non-exotic spec¬ 
tra. Gauge field dynamics are modeled via flux tube exchange factors. The energy determined for 
the ground states with J K = ( + ) is 2.22 GeV (2.50 GeV). A lower energy negative parity state 

is consistent with recent lattice results. 


A system of interacting, non-relativistic constituent quarks is the most simple, realistic 
model of hadronic systems. Solving the many-body Schrodinger equation to determine 
wave functions within this simple model is still a formidable task owing to the strong 
flavor, spin, and color dependence of the quark-quark interaction and traditionally re¬ 
quires some array of approximate methods to solve it. The controversial status of the 
recent evidence of a flavor exotic five-quark state warrants a careful treatment of this 
strongly interacting system. We have determined the wave function of five interacting 
constituent quarks in the flavor-exotic multiquark hadronic sector is solved using the 
variational Monte Carlo (VMC) technique |l]. This technique is known to yield upper 
bounds on the ground state energy accurate to the level of a few percent in light nuclei 
with the number of nucleons A < 6 ||2] • 

Recent experimental evidence (for reviews see [3;, 11) has revived interest in the mul¬ 
tiquark flavor-exotic sector of the hadronic spectrum. Various model calculations have 
been reported (an overview is given in |5]) which study the existence of a strangeness 
S — +1 resonance, dubbed 0 + , about 100 MeV above threshold to nucleon-kaon decay 
with low mass of 1540(2) MeV and possibly extremely narrow width of 0.9(3) MeV 
[6]. Lattice results are available, including Refs. fli . High statistics photoproduction 
data of the reaction yp —> K°K + n from the CLAS collaboration H] sets an upper limit 
on the yield of the 0 + relative to the A*(1520) yield at 0.2%. No definitive structure in 
the nK -invariant mass spectrum is observed at 1540 MeV in this experiment. No experi¬ 
mental information is available on the spin or parity of the state in any experiment done 
to date. 

We study the 9 + in the non-relativistic flux tube quark model with one-gluon ex¬ 
change (OGE) and one-pion exchange (OPE) (between the light quarks). We take se¬ 
riously the possibility that 0 + is a narrow resonance and calculate its mass as a stable 


TABLE 1 . Variational energies of the states of pentaquark 
states in GeV. The mass is reported as M g +(n;J%) = (T) + 
(V) — 385.5(5) MeV. Statistical errors for M d + are < 5 MeV. 


n\J% 

l;H 

1;1 + 

2;1+ 

3;0+ 

3; 1 + 

4;0+ 

4;1 + 

M d + 

2.22 

2.50 

2.57 

2.75 

2.81 

2.83 

2.88 

(f) 

1.68 

2.13 

2.02 

2.03 

2.00 

1.92 

1.90 

<v> 

0.92 

0.74 

0.93 

1.10 

1.19 

1.29 

1.36 


state with respect to strong interactions. The approach adopted here will be to work with 
a general, completely antisymmetric wave function and determine dynamically which 
flavor-spin-color-orbital (TSCL) structures are favored in the constituent quark model 
(CQM). The model Hamiltonian used in this work is fixed by the single hadron spec¬ 
trum Em an d the six-quark (deuteron) properties Eh. It is given as a sum of kinetic 
energy, pair potentials for OGE and OPE, and confining terms. The two-body operator 
potentials for OGE (OPE) interactions are determined in the non-relativistic reduction 
of the tree level amplitude for qq or qq scattering. We use a many-body confining inter¬ 
action based on the “flux tube” model Eh. We consider flux tube topologies (there are 
three) for the 4qq states consistent with gauge invariance. Matrix elements of operators 
which connect different flux tube topologies are modeled by a quark position dependent 
factor which models the dynamics of the gauge fields Eh. 

The state vector for the five quarks, including two-body isospin, spin, color and spatial 
correlations is given by 

(d 

i<j 

where the symmetrized product of two-body correlation operators Fy acts on an uncor¬ 
related TSCL state, The state |4>s) and, therefore, is completely antisymmet¬ 
ric under exchange of the coordinates of the u and d quarks. 

We considered a large set of states obtained by considering states of total isospin 
T — 0, S — 0,1, and L = 0.1 to yield states of negative and positive parity. The negative 
parity state |1;1“) corresponds to the state of 4 q with positive parity, 714 — +1 with 
zero orbital angular momentum. The positive parity states, \n = 1,... ,4;/^) correspond 
to states of 4 q with negative parity and include one unit of orbital angular momentum. 
Here J% denotes the spin of the four light quarks, J 4 and the parity, K of the pentaquark 
state. 

The evaluation of expectation values of the variational wave function I'Ps) are effected 
via Monte Carlo integration of the 15 dimensional space R = (it, ... ,r 5 ). Results for the 
variational energies of negative and positive parity states are shown in TableU] The mass 
of the state is reported as 


M 0 + (n;J%) = 4 m q +m s + (t) + (V) - V 0 (4qq) (2) 

where ( T ) and (V) are the variational kinetic and potential energies, respectively, for the 
states considered. The quark masses are taken as m q = 313 MeV, m s = 550 MeV and 
Vo is a constant fit to the nucleon mass. In reality, the true positive parity ground state 
is some admixture of the six positive parity states. However, the lowest energy state 





dominates since the states \n\J%) are orthonormal: applying the correlations as in Eq.© 
introduces off-diagonal elements which are quite small. States with n ^ 1 contribute 
< 10% to the true positive parity ground state. 

Comparison of the CQM with lattice results at large pion mass is meaningful since the 
effects of chiral symmetry are negligible. We note that our negative parity state 11; 1 _ ) 
lies below the positive parity state 11; 1 + ), consistent with the lattice data in Refs. O, ISO- 
This is true despite the fact that the positive parity state has much stronger attraction 
from OGE and OPE contributions, attributable to the high degree of symmetry of the 
TS state (see Table©. The (V) of the negative parity state is 180 MeV below the lowest 
positive parity state. Exciting one unit of orbital angular momentum on the other hand 
raises the kinetic energy significantly, about 450 MeV. 

We have shown that within the non-relativistic flux tube CQM with OGE and OPE 
interactions the negative parity state is lower than the positive parity state by ~ 280 
MeV when the 0 + is assumed to be a narrow resonance. Though the higher symmetry 
of the lowest lying positive parity state significantly decrease its potential energy, its 
unit excitation of orbital angular momentum raises the energy above that of the negative 
parity state. Quark correlations are crucial in making these determinations and should 
not be ignored. 
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